Formation of soliton trains in Bose-Einstein condensates by temporal Talbot effect 
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We study the recent observation of formation of matter-wave soliton trains in Bose-Einstein 
condensates. We emphasize the role of the box-like confinement of the Bose-Einstein condensate 
and find that there exist time intervals for the opening the box that support the generation of 
real solitons. When the box-like potential is switched off outside the existing time windows, the 
number of peaks in a train changes resembling missing solitons observed in the experiment. Our 
findings indicate that a new way of generating soliton trains in condensates through the temporal, 
matter-wave Talbot effect is possible. 
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The notion of soliton belongs to the most popular con- 
cepts in physics. Any localized and long-living structure 
is readily called a soliton. Solitons are identified with so- 
lutions of nonlinear equations and appear in many con- 
texts of science ranging from physics (nonlinear optics, 
hydrodynamics, particle physics, etc.) to molecular bi- 
ology. Solitons are formed because of the existence of 
nonlinear interaction in the system which cancels the dis- 
persion and hence allows for the propagation of shape 
preserving objects. 

In the case of dilute atomic quantum gases the nonlin- 
earity is determined by the effective interaction between 
atoms that can be both repulsive or attractive. For repul- 
sive uniform condensates the appropriate nonlinear equa- 
tion (i.e. the Gross-Pitaevskii equation) predicts dark 
soliton (a hole in the density associated with a phase 
jump) as a solution Jj. Such excitations of Bosc-Einstcin 
condensates have been already observed in experiments 
with trapped alkali atoms Pj . 

Generating bright solitons in atomic quantum gases is 
more difficult task because it requires working with at- 
tractive condensates. Due to a collapse, in such samples 
the number of atoms is limited and small (see Ref. |3]). 
This obstacle has been overcome in two ways. In the first 
one, a large repulsive condensate is formed and then the 
interactions are changed from repulsive to attractive by 
using the Feshbach resonance technique 0, 0| , whereas in 
the second attempt j^d the optical lattice and the notion 
of negative effective mass are utilized. 

The goal of this Letter is twofold. First, we focus on 
the experiment of Ref. 4] where the trains of bright 
solitons are generated in attractive condensate of ^Li 
atoms. Existing theoretical explanations involve quan- 
tum phase fluctuations :7] or modulational instability 
1^ as the main mechanism leading to the observed 
structures. None of these papers, however, attempt to 
model the experiment of Ref. |j| . In our Letter we report 
on calculations that correspond to the same values of all 
parameters as in Ref. Q with realistic three-dimensional 
geometry. Secondly, we concentrate on the role of the 
box-like confinement and show that there exists a regime 
where multiple scattering of gas from the walls of the box 
might result in formation of "real" solitons, i.e. multi- 
peak structures that undergo elastic collisions. We sug- 



gest a new version of the experiment and specify physical 
conditions supporting the formation of "genuine" soliton 
trains. The underlying physics of our proposition can be 
understood with the help of the temporal Talbot effect, 
already experimentally observed in the context of Bose- 
Einstein condensate of sodium atoms diffracted by the 
pulsed grating [lOj. 

To describe the experiment of Ref. Q we start, as 
the Ref. [1| does, with the time-dependent dissipative 
Gross-Pitaevskii equation 
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where 'ip{r,t) is the macroscopic wave function of Bose- 
Einstein condensate of ""Li atoms, Vtr{z, p) — m {cu'^z'^ -\- 
/2 + Vbox is the axially symmetric trapping poten- 
tial, g = ^Trh^as/m (with as being the scattering length 
which determines the strength of the interaction). The 
number of atoms equals iV = 3 x 10^, the geometry of 
the harmonic confinement is given by Uz = 27r x 4 Hz 
and ijo± = 27r X 800 Hz, and the box-like potential Vbox 
(the "end caps") of length of L = 12 osc. units is po- 
sitioned on the slope of the harmonic potential, shifted 
by 15 osc. units from its center. The initial value of the 
scattering length of ^Li condensate prepared in the end 
caps, which is positive and equals 200 ao (ao - Bohr ra- 
dius), is changed within approximately 10 ms to its final 
value Os — — 3ao and then the system is kept in the box- 
like potential for further 17ms. Finally, the end caps are 
turned off. As opposed to the Ref. our calculational 
parameters are the same as in the experiment of Ref. . 

The imaginary term in Eq. describes the losses 
due to three-body recombination processes ^J, 0. Its 
presence is necessary to get the agreement with exper- 
imental results. Since the losses were not investigated 
experimentally, we follow the references just mentioned 
and put 7 = 2.05 x 10^^^ cm^ s~^. The dissipative term is 
turned on, according to the observation in 4], when the 
interaction strength becomes negative and is strongly di- 
minished (by a factor of 100) when the end caps are off. 
When the losses are kept constant during the whole cal- 
culations the final number of atoms is much less than the 
number observed in the experiment. On the other hand. 
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FIG. 1: Moving soliton train of ^Li atoms near the turning points and at the center of the oscillation. Successive images 
correspond to 6 ms, 48 ms, 125 ms, and 185 ms after the end caps are off. The number of atoms, the trap parameters, and the 
way the scattering length is changed are the same as in the experiment of Ref. jj. Color images are available at [T^. 



decreasing the dissipative term leads to the condensate 
collapse. Another way of getting agreement with experi- 
mental data (i.e. the final number of atoms) is to lower 
the initial number of atoms while keeping the same value 
of 7 all the time. Hence, the dissipative term plays a 
crucial role in stabilizing the system. Certainly, it has an 
influence on the collapse studied recently in Ref. 9]. 

After switching off the box potential the bosonic cloud 
starts to oscillate in harmonic trap. As it is illustrated 
in Fig. n the clouds actually breaks into several peaks 
which propagate in the potential for many oscillatory cy- 
cles. However, we observe that certain peaks disappear 
during the evolution and reappear again later. For exam- 
ple, in the second frame of Fig. ^when the condensate 
goes through the oscillatory center from the left to the 
right, the number of distinguishable peaks is 3, whereas 
in the lowest frame (the system goes now to the left) that 
number is increased to 4. Such structures can not be, in 
fact, considered as solitons. 

Explaining the experimental results of Ref. we, as 
opposed to other theoretical works 0j H, 01 j emphasize 
the role of the box-like potential. We demonstrate how 
important is the time when it is open and its location. 
To this end, we have performed calculations for a smaller 
sample of = 10^ ^Li atoms created with initial scat- 
tering length Us ~ lOOao, which is next changed within 
10 ms to the final value = — 3ao. The dissipative term 
is kept constant all the time. Our results are summarized 
in Fig. 121 The end caps are positioned symmetrically 
with respect to the center of the harmonic trap. Only 
in such a configuration we discover the existence of time 
windows, i.e. the time intervals for the opening the end 
caps that support the generation of train of real solitons 
(multipeak structures with preserved number of peaks). 
This is illustrated in Fig. |21 where the size of the box- 
like potential is L = 4.0 osc. units (75.6 /im). Frame 
(a) shows the density in the case when the end caps are 
switched off within the time window (61 ms after the in- 
teraction strength is changed) . Here the number of peaks 
is still the same. On the contrary, when the end caps are 
off at the time outside the time window, the number of 
peaks changes resembling the missing solitons observed 
in the experiment 0|. Calculations also show the exis- 
tence of further time windows, especially for shorter end 
caps. The duration of the time window is about 3 ms and 
increases when the size of the end caps gets larger. 



It is important to understand the role of the location of 
the box-like potential. We have checked numerically that 
the time windows persist even if the end caps are shifted 
from the center of the harmonic potential by several per- 
cent, depending on the size of the end caps. Moreover, we 
do not see the time windows when the size of the end caps 
is too large. It is clear that the discussed phenomenon is 
somehow spoiled by the presence of the harmonic trap. 
As will be explained later, the disappearance of time win- 
dows is caused by the loss of Talbot-type recurrence due 
to nonuniform shift of energy levels of the box poten- 
tial forced by the harmonic trap. Certainly, the most 
favorable conditions for the existence of time windows 
are those when there is no axial harmonic confinement. 
Such calculations have been already performed, see Fig. 
1 in Ref. Q. However, the authors of Ref. Q have not 
investigated carefully the influence of the delay time be- 
tween the removal of the end caps and the switching off 
the scattering length and they overlooked the possibility 
of existence of time windows. In Fig. |3]we plot the axial 
density profiles in the case when the box-like potential is 
turned off just at the time when the scattering length is 
changed (left column) and in the case when the box is 
removed within the time window. The left column shows 
that, in fact, the condensate is split into 6 peaks (not 5 
as reported in Q) and its evolution is violent in a sense 
that the number of peaks changes for a long duration. On 
the other hand, as demonstrated in the right column, the 
number of peaks (solitons) is settled much earlier when 
the end caps are switched off within the time window. 

The origin of time windows is due to the temporal 
Talbot effect. Therefore, we start with a purely linear 
case and consider the evolution of the symmetric wave 
packet in an "open" rectangular one-dimensional box. 
The walls of the box are very high and therefore we will 
expand the initial wave packet in the basis of the symmet- 
ric eigenstates of infinitely deep well potential that are 
\/2/L cos {knx) with fc„ = r(2" + l)i where n = 0, 1, ... 
and L is the length of the well. Note that the wave func- 
tions are zero in the classically forbidden region. The 
time evolution reads 

V'(x,t) = ^a„e-'^''(^"+i)'*/''y|cosfc„x, (2) 

n 

where Eq = fi.^7r^/(2mL^) and {a„} are the coefficients 
of the expansion of ^{x,0). Fourier transform of |2Jl is 
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FIG. 2: Illustration of the importance of a time when the end 
caps are switched off. Frame (a) shows the condensate density 
at 30 ms, 66 ms, and 108 ms after the box is off, i.e. 61 ms 
after the interaction strength is changed. This is the case of 
the time window when, as a result, the structure consisted of 
the same number of peaks is formed. As opposed, frame (b) 
resembles Fig. Here, the number of peaks changes during 
the evolution. Successive snapshots correspond to 30 ms, 66 
ms, and 139 ms after the end caps are turned off outside 
the time window (52 ms after the change of the interaction 
strength). 



given by 
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According to the above formula the momentum distribu- 
tion is recovered after a period of Trev = 7r/4 (in units of 
h/Ef)) and this is the result of the Talbot-type recurrence 
occurring in the phase factor ( "n(n-\- 1)" dependence on 
the quantum number n). However, this formula exhibits 
more structure. At time T^"n = tt/S the resulting mo- 
mentum distribution forms a set of fully separated groups 
centered at momenta k = kn- This can be verified (e.g. 
numerically) taking into account several facts: (1) the 
form of the phase factor in Eq. 0; (2) the localization 
of the function sin(x)/x around x = 0; (3) assumed weak 
dependence of coefficients a„ on n. 

If the box-like potential is switched off at times around 
'^win — ""/^ groups of momenta transform into sepa- 
rated wave packets which just spread out during the fur- 
ther evolution. This spreading can be stopped by turn- 
ing on the nonlinear term and then the wave packets are 
transformed into Zakharov solitons p|. Hence, the non- 
linearity is essential to built the solitons from the sep- 
arated wave packets. Here, the separated wave packets 
originate from the temporal Talbot effect combined with 
the appropriate initial conditions. In a real experiment 
initial conditions are formed as a result of modulational 
instability driven by the condensate collapse. 
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FIG. 3: Axial density profiles of a condensate of 10* ^Li atoms 
during the evolution under no axial confinement. The size of 
the end caps equals 189 /im. The left column corresponds to 
Fig. 1 in Ref. Q when time window condition is not fulfilled. 
For the right column the end caps are turned off within the 
time window (314ms after the scattering length is changed). 
In both cases the time t = means the time when the end 
caps are removed. 



Looking at the formula |j3Jl, it is clear that the presence 
of the box-like potential is essential since it allows the 
condensate wave function for multiple reflexions from the 
walls resulting in the interference pattern found above. 
However, "small addition" of the harmonic potential does 
not spoil the picture. This can be proved by the pertur- 
bation calculus. It turns out that for not too large size of 
the box L (L < 5 osc. units) all the levels but the lowest 
one are shifted approximately uniformly. In such a way 
the new time scale that is proportional to the recipro- 
cal of the difference of the shifts of the two lowest levels 
(and approximately equals the trap period) can be intro- 
duced. Based on the formula Q only those time windows 
survive that occur at times shorter than the trap period 
(although "revivals" are observed for long enough times). 

We can also prove that the temporal Talbot effect sur- 
vives under the presence of the weak nonlinearity. For 
small enough attraction the nonlinear term that appears 
in the Gross-Pitaevskii equation can be treated as a per- 
turbation to the linear case. It turns out that all single- 
particle levels are shifted up uniformly. Therefore, from 
the formula ©, one should expect formation of groups 
in momentum space at the same times and of the same 
duration as in the linear case. 

The temporal Talbot effect persists also when the non- 
linearity is larger and cornparable to that one present in 
the experiment of Ref. [J. As it is the case, we plot 
in Fig. 0] the axial momentum densities (at the time 
when the end caps are off) corresponding to the frames 
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FIG. 4: Spatial (left column) and momentum (right column) 
densities at the time when the end caps are off. The upper 
frames illustrate the case of the time window in which the 
momentum density shows well developed groups of momenta. 
These groups are next transformed to groups in a position 
space. The lower frames correspond to the case when there are 
no separated peaks in a momentum space and consequently 
the probability flow leads to "missing solitons" structures. 

(a) and (b) of Fig. |2 Fig. 0] proves that the time win- 
dows survive under the conditions when both the box- 
like and harmonic confinements and the nonlinearity are 
present. When the end caps are off just within the time 
window (the upper frames) the momentum density shows 
five distinguishable peaks whereas the density in the po- 
sition space consists of two broad peaks with partially 
developed three subpeaks. Upper frames in Fig. |2] con- 
firm that afterward five (not two nor six) solitons are 
developed. It means that all momenta groups have been 
transformed to peaks in position space. These peaks os- 
cillate with the trap period, collide when they meet at 
the center of the trap and then reappear. Their motion 
is particle-like. If one considers a point-like particle mov- 
ing according to the Newton equation, initially placed at 
the trap center and assign the initial velocity determined 
by the maximum value of the momenta group, the par- 
ticle will follow the density peak. Hence, it is reasonable 
to use the name solitons for the density peaks. 

When the end caps are switched off outside the time 
window (lower frames in Fig. 0)) the number of density 
peaks changes during the evolution (see the lower part 



of Fig. 12). This is because the momentum density does 
not consist of well separated groups and the probability 
can flow from one group to the other forming the "missing 
solitons" structures. This is, in fact, the case of numerical 
calculations reported in Ref. 0. Fig. 5 of Ref. Q 
shows that the number of peaks changes in time while 
the condensate is axially confined. On the other hand, 
when the axial confinement is off the condensate always 
ends in a state with a fixed number of solitons. 

The time the first time window appears is slightly 
larger than T^™j — tt/S predicted from the linear theory. 
We found that it is shifted by an amount of the order of 
characteristic time scale H/ (gN/V), which in our case is 
equal to several milliseconds. 

In conclusion, we have shown the importance of the 
box- like potential in experiments like that of Ref. 
Contrary to missing solitons structures genuine bright 
solitons can be generated in a repulsive condensate by 
changing the sign of the scattering length with the help 
of the Feshbach resonance technique and keeping the con- 
densate in the end caps for long enough time. In fact, 
what we propose to create solitons is to utilize the Tal- 
bot effect that is well known from the linear physics. It 
turns out, however, that this effect survives under the 
condition of weak nonlinearity what was proved by us- 
ing the perturbative analysis but also under the condi- 
tions of the experiment of Ref. as was demonstrated 
by direct solving of the Gross-Pitaevskii equation. Per- 
haps, it is not surprising since the temporal Talbot effect 
was already experimentally observed for a condensate of 
sodium atoms, although in a different situation when the 
condensate was diffracted by a pair of pulsed gratings and 
the periodicity with respect to the time delay between 
pulses in a series of condensate images was discovered 
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